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ABSTRACT 

Aims. Tau is a detached and single-lined interferometric-spectroscopic binary as well as the most massive binary system of the 
Hyades cluster The system revolves in an eccentric orbit with a periodicity of 140.7 days. Its light curve furthermore shows a complex 
pattern of 6 Scuti-type pulsations. The secondary has a similar temperature but is less evolved and fainter than the primary. In addition, 
it is rotating more rapidly. Since the composite spectra are heavily blended, the direct extraction of radial velocities over the orbit of 
component B was hitherto unsuccessful. Our aim is to reveal the spectrum of the fainter component and its corresponding Doppler 
shifts in order to improve the accuracy of the physical properties of this important "calibrator" system. 

Methods. Using high-resolution spectroscopic data recently obtained with the Elodie (Observatoire de Haute-Provence, France) 
and Hermes (Roque de Los Muchachos, La Palma, Spain) spectrographs, and applying a spectra disentangling algorithm to three 
independent data sets including CfA spectra (Oak Ridge Observatory, USA), we derived an improved spectroscopic orbit. We next 
used a code based on simulated annealing and general least-squares minimization to refine the orbital solution by performing a 
combined astrometric-spectroscopic analysis based on the new spectroscopy and the long-baseline data from the Mark III optical 
interferometer 

Results. As a result of the performed disentangling, and notwithstanding the high degree of blending, the velocity amplitude of the 
fainter component is obtained in a direct and objective way. Major progress based on this new determination includes an improved 
computation of the orbital parallax (still consistent with previous values). Our mass ratio is in good agreement with the older estimates 
of Peterson et al. (1991, 1993), but the mass of the primary is 15-25% higher than the more recent estimates by Torres et al. (1997) 
and Armstrong et al. (2006). 

Conclusions. The evolutionary status of both components is re-evaluated in the light of the revisited properties of 9^ Tau AB. Due to 
the strategic position of the components in the turnoff region of the cluster, the new determinations imply stricter constraints for the 
age and the metallicity of the Hyades cluster We conclude that the location of component B can be explained by current evolutionary 
models, but the location (and the status) of the more evolved component A is not trivially explained and requires a detailed abundance 
analysis of its disentangled spectrum. The improved accuracy (at the 2% level) on the stellar masses provides a useful basis for the 
comparison of the observed pulsation frequencies with suitable theoretical models. 

Key words, astrometry - techniques: high angular resolution - stars: binaries: visual - stars: binaries: spectroscopic - stars: funda- 
mental parameters - stars: individual: 6- Tau 



1. Introduction 

A modern research topic that we are cuiTently pursuing is the 
study of binary and multiple stars with at least one pulsating 
component. The advantages of studying pulsating components 
in well-detached systems are manyfold, e.g. both the theories of 
stellar evolution and of stellar pulsation can be accurately tested 
and refined. Accurate component properties compared to suit- 
ably chosen theoretical isochrones indeed allow to obtain infor- 
mation on the object's age and evolutionary path and are usually 
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* Based on observations obtained at the 1.93-m telescope of the 
Obserx'atoire de Haute Provence, the 1.2-m Mercator telescope at the 
Roque de los Muchachos Observatory (in the framework of the Hermes 
Consortium) and the 1.5-m Wyeth telescope at Oak Ridge Observatory. 



necessary to help discriminate among various possible pulsation 
models. 

We selected the 6 Scuti star 0^ Tau for a careful study for 
three reasons: (a) it is a detached, spectroscopic binary re solved 
by long-baseline interferometry (lArmstrong et al.l [20061 here- 
after AM06); (b) it is a member of the Hyade s open cluster at a 
mean distance of 45 pc ("Perrvma n et al.|[T998l) : (c) the evolution- 
ary status of the pr imary component is still under heavy de bate 
(lTorresetal.1119971 hereafter TSL97; iLastennet et aPl 19991 and 
AM06). 

6)2 Tau (HD 28319 = HIP 20894 = 78 Tau) is the most mas- 
sive main-sequence star of the Hyades cluster: it is located at the 
main-sequence turnoff region of the isochrone best fitting the in- 
dividual members of the cluster (Lebreton et al. 2001, hereafter 
LEOl). It forms a quadruple system with 0' Tau (= 77 Tau), a 
common proper motion companion at an angular separation of 



1 



Torres, K. et al.: Orbit, orbital parallax and masses of Tau AB 



5. '6 . The brighter component is a single-lined spectroscopic bi- 
nary (SBl) which was resolved by long-baseline interferometry 
(Shao et al. 1991; Pan et al. 1992; Hummel & Armstrong 1992). 
The orbital period of the binary is 140.7 days with an eccentricity 
of about 0.7 . The primary component, 0^ Tau A, has been clas- 
sified as A7 III and rotates with Vsin / = 70 km s"^ dFremat et alj 
12006-^. It is very hard to detect the secondary, 9^ Tau B, spec- 
troscopically because the Doppler shifts are only a fraction of 
the width of its broad spectral lines. The secondary is less 
evolved and therefore fainter than 6^ Tau A. TSL97 treated 
9^ Tau as a double-lined spectroscopic binary (SB2) using a 
2D-cross correlation method in which they considered the extra 
pull of the secondary in order to obtain improved radial veloc- 
ities of the primary component. However, they were unable to 
obtain reliable radial velocity measurements for 9^ Tau B and 
they only observed its (orbital) influence on the velocities of 
the brighter companion. For these reasons, also the mass ratio 
could not be directly determined. Still, they were able to de- 
rive a interferometric-spectroscopic orbit and to determine the 
component masses and the distance of the binary by exploring 
and assuming a range of values for the mass ratio and the rota- 
tional velocity of the secondary star. The orbital parallax of their 
solution agrees well with the Hipparcos trigonometric parallax 
(lPerrvman& ESAl[T997l) . The outcome is that the components 
have different projected rotational velocities: TSL97 obtained a 
best fit assuming = 38 km s"' and VfiSin/ =110 km s"', 
with resulting component masses of - 2.42 + 0.30 M© and 
Mb = 2.11 +0.17 Mq. Recently, the component masses and lu- 
minosities have been redetermined by AM06. From their inter- 
ferometric data coupled to the Hipparcos secular (i.e. proper- 
motion based) parallax (de Bruiine et al. 2001) and a compati- 
ble choice of spectroscopic orbits, these authors obtained com- 
ponent masses of Ma = 2.15 + 0.12 Mq and Mg = 1.87 ± 0.11 
Mo. 

As a member of a well-studied cluster, both the metallicity 
and the distance of 0^ Tau are known within narrow boundaries. 
Such a large amount of information concerning the fundamental 
properties of the components of a detached binary system (with 
components in a different evolutionary phase) has led to vari- 
ous attempts of c onfrontation with stellar evolutionary models . 
Before Hipparcos dPerrvman & ESAll997h . lKr61ikowskal(ll992h 
concluded that the evolutionary status of Tau was either in 
the thick H(ydrogen)-shell burning phase (without overshoot- 
ing) or in the H-core burning phase (with overshooting). From 
the location in a colour-magnitude diagram and a best fit with an 
isochr one of age ^ 630 My r and metallicity Z = 0.027, TSL97 
(using iBertelli et al.l d 19941) 's models) concluded that the pri- 
mary is in a phase near H-core exhaustion, immediately preced- 
ing the phase of overall contraction. But, because both binarity 
and fast r otation may afi'ect the colour indices, the uncertainty 
remained. lLastennet et aljfl999l) used the binary to evaluate var- 
ious stellar evolution models stating that "the three theoretical 
models allow to fit correctly the system (...) in agreement with 
the more recent constraints available about the metallicity of the 
Hyades cluster" . They found that the primary could be in the H- 
core burning (end) phase or in the H-shell burning (beginning) 
phase depending on its metal lici t y. In a n extensive study of the 
Hyades cluster, Ide Bruiine et alJ (l200ll) again used 9^ Tau and 
concluded that the agreement with the CESAM isochrones (in- 
cluding convective core overshoot) was remarkably good. LEOl 
derived a maximum age of 650 Myr and the initial helium con- 
tent of the cluster by comparing the mass-luminosity relation 
based on a set of five binaries also containing 0^ Tau with predic- 
tions from models appropriate for the Hyades (using Y = 0.26). 



They issued a warning for the stars located in the turnoff region 
claiming that "the interpretation is complicated by the effects of 
rotation and overshooting that make either model or photometric 
data uncertain" and that "improvement of the mass of 0^ Tau A 
would certainly better constrain the overshooting by anchor- 
ing the star more precisely on its isochrone" . More accurate 
masses is precisely what AM06 obtained. However, the lower 
masses and/or brighter luminosities do not c onform with some of 
the recent stellar evolutionary models, i.e. iGirardi et alj d200d) 
and LEOl, f or the age and metallicity of the Hyades. Finally, 
Yild iz et al.l (|2006) tried to model both components of 0^ Tau us- 
ing models including differential rotation. They concluded that 
good agreement was found for 9^ Tau A with Z = 0.024 and an 
age of 700 Myr However, the same model could not explain the 
location of the secondary component in the colour-magnitude di- 
agram. The evolutionary status of both components of 0^ Tau is 
therefore still a non-trivial issue. 

The system has one more attractive feature; its primary com- 
ponent is a typical multiperiodic 6 S cuti star. Var ious multi- 
site campaigns have been conducted. iBreger et al.l (1198 9) ob- 
tained five closely spaced and stable frequencies, all of which 
had amplitudes below 0.01 mag. They discarded rotational split- 
ting since it could not explain the observed frequency separa- 
tion s and proposed a mixt ure of modes of different I and m val- 
ues. Kennellvet all d 199 6*) discussed a large set of radial velocity 
and line profile data from which up to seven frequencies emerged 
with only three frequencies in common with the previous analy- 
sis. They suggested long-term (> 6 yr) amplitude variability and 
a combination of low and high degree modes. Amplitude vari - 
ability on a 10 yr time scale is also claimed by Li et al. (1997|). 
Both components are located within the lower Cepheid insta- 
bility strip (where the 6 Scuti stars are observed) but it seems 
well established that the more massive prima ry is the pulsating 
star (Breger et al. 1989; Kennelly et aLlll996h . The results of a 



vast multi-site campaign (IBreger et a 120021) revealed the pres- 
ence of 1 1 frequencies in the ran ge 10-15 c/d and 2 frequencies 
in the range 26-27 c/d (see also lPoretti et al.ll2"002l) . After hav- 
ing modelled the light-time corrections, they attributed the first 
1 1 frequencies to pulsations of the primary component while the 
last two frequencies probably originate from the secondary com- 
ponent. Currently, because of the limited frequency range of the 
detected modes, the complex frequency spectrum does not allow 
a proper mode identification. 

Because they are located in the turnofi" region of the clus- 
ter, the exact locations of both stars in the H-R diagram allow 
to distinguish between different evolutionary models and differ- 
ent isochrones, thus making them extremely useful for the con- 
straints they impose on the chemical composition and the age of 
the Hyades (provided that their physical properties can be de- 
termined very accurately). As "calibrator stars", they could in- 
dicate how important the mixing processes (such as convective 
core overshooting or rotational mixing) are in the internal layers 
of stars from the Hyades cluster (LEOl). The knowledge of accu- 
rate fundamental component properties furthermore holds good 
potential for a reliable pulsation modelling of the two stars. 

The outline of this article is the following: the observations 
are presented in Sect. 2 while the analysis method of spectra 
disentangling is explained in some detail in Sect. 3. In Sect. 4, 
we will compare the component spectra obtained from the ap- 
plication of the spectra disentangling method with an observed 
spectrum of a reference star as well as with synthetic spec- 
tra. Especially in the case of the fainter component, not easily 
discernable in the observed composite spectra, this will show 
that the mathematically reconstructed spectra are plausible from 
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the viewpoint of physics. Sect. 5 deals with the orbital analysis 
obtained by combining both spectroscopic and interferometric 
data. In Sect. 6, we thoroughly discuss the quality of our results 
and their implication in the light of the component's evolution- 
ary stages. We end by mentioning the new perspectives and the 
planned future work. 

2. Observations 

The spectroscopic data consist of (a) 44 echelle spectra 
(R=42000; S/I^CK200-300 in V; resolution element of 7.1 
km s ') obtained by us with Elodie at the 1.93 m telescope of 
the Observatoire de Haute-Provence covering the wavelength 
range from 389.5 to 681.5 nm; (b) 13 echelle spectra (R~85000; 
S/N~220-340 in V; resolution element of 3.5 km s"') obtained 
by us with the Hermes spectrograph at the 1.2 m Mercator tele- 
scope located at the Roque de Los Muchachos Observatory, La 
Palma, covering the wavelength range from 377 to 900 nm; (c) 
70 echelle spectra (R~35000; S/N~50; resolution element of 8.5 
km s~') obtained by TSL97 using the Center for Astrophysics 
(Cf A) spectrograph mounted on the 1.5 m Wyeth reflector at the 
Oak Ridge Observatory, covering a wavelength range of 26 A 
around 519 nm; (d) 16 single-order coude spectra (R~25000; 
S/N~100-200; resolution element of 14 km s"') obtained at the 
2-m telescope of the Ondfejov Observatory of the Astronomical 
Institute of the Academy of Sciences of the Czech Republic, cov- 
ering the wavelength range 517-589 nm. 

The Elodie data w ere reduced using the Intertacos pipeline 
(iBaranne et al.l Il996h available at the telescope, while the 
Hermes data were treated using the Hermes reduction pipeline 
(iRaskin et al.ll20ld l). These reduction procedures perform the or- 
der extraction, the offset and flat-field correction, and the wave- 
length calibrations. The resulting wavelength scale was cor- 
rected for the Earth's motion relatively to the barycenter of the 
solar system using the Iraf software package. The spectra col- 
lected at the Ondfejov Observatory were reduced with the pro- 
cedure described by , Skoda & Slechta (2002). 

In total, 127 spectr^ covering the entire orbital cycle and full 
range in radial velocity amplitude were available for the anal- 
ysio We furthermore made use of 34 best-fit angular separa- 
tions (p) and position angles (0 ) derived from the measurements 
obtained by [Armstrong et alj (Il992h . iHummel & ArmstrongI 
(119921) and AM06 with the Mark III long-baseUne optical in- 
terferometer. 

3. Spectra disentangling 

We used the spectra dise ntangling code F DBinary V.3 
(Rel. 30.01.09) developed bv llhiic et all (|2004 . which deter- 
mines the individual contributions of the components to the com- 
posite spectra together with the orbital parameters in a self- 
consistent way. The FDBinary runs are based on the input of 

' From hereon, S/N = S/N per resolution element 

^ The 44 spectra from the Elodie spectrograph can be retrieved di- 
rectly from the Elodie database (http://atlas.obs-hp.fr/elodie/ 1. The 70 
spectra from the Harvard Smithsonian Center for Astrophysics (CfA) 
and the 13 ones from the Hermes spectrograph are available upon re- 
quest from the authors of the paper Torres et al. 1997 and from us, 
respectively. 

^ Unfortunately, the Ondfejov single-order coude spectra did not en- 
ter the subsequent disentangling procedure due to their lower resolution 
and lower S/N (caused by unfavourable weather conditions during the 
observations) compared to the sets of echelle spectra. 



Table 1. Spectral intervals used in the spectra disentangling anal- 
ysis. 



Index 


4' 

(A) 


\f 
(A) 


N 




Rl 


4660.37 


4725.97 


2363 


0.364 


R2 


4759.53 


4961.50 


4520 


0.364 


R3 


5178.81 


5292.47 


3700 


0.360 


R4 


5421.21 


5475.42 


1700 


0.359 


R5 


5520.36 


5610.08 


2700 


0.359 


R6 


6090.67 


6127.51 


1035 


0.355 


R7 


6127.55 


6187.31 


1665 


0.353 



Notes. Aj and Af refers to the initial and final wavelength in Angstroms; 
N refers to the number of bins contained in the observed composite 
spectra; and {bI^a corresponds to the light ratio of the components de- 
termined as explained in Sect.|3] 



observed spectra at epochs which, in an ideal case, uniformly 
cover the orbital radial velocity range. The observed spectrum is 
assumed to be the combination of two time-independent compo- 
nent spectra that are shifted in wavelength with respect to each 
other (according to the derived Keplerian orbit). The code uses 
the multi-dimension al, non-linear opti mization technique of the 
(downhill-)simplex dPress et al.lll992l) to determine the orbital 
parameters, whereas the intrinsic component spectra are com- 
puted by the algorithm of spectra separation (using the radial ve- 
locities estimated from Kepler's equations). Commonly, many 
runs starting from various points in a subspace of the orbital pa- 
rameter space are launched, and different sizes for the initial sim- 
plex can be explored. Convergence is attained when the size of 
the simplex shrinks below a specified level in a specified max- 
imum number of iterations. The separation algorithm works on 
the Fourier components of the observed spectra using singular 
value decomposition for quasi-singular sets of equations. 

Our analysis was divided in two steps: (a) the determination 
of the orbital parameters from a selected wavelength region and 
(b) the computation of the intrinsic component spectra, keeping 
the orbital parameters fixed during the convergence, for other 
wavelength regions (cf. Table [T]|. 

In step (a), the search for the orbital parameters was done 
using 127 spectra weighted according to their S/N ratio in the 
wavelength region R3 (cf. Table [TJ. This region was chosen 
for it allows the highest accuracy on the determination of the 
component radial velocities among all selected regions (cf. 
Fig. [1] the standard deviation shows a deep local minimum in 
that region) . Since the CfA spectra cover the large range of 
radial velocity amplitude near periastron passage but have a 
smaller wavelength range of 26 A, we only used a part of R3. 
Moreover, the CfA, Elodie and Hermes spectra together provide 
sufficient spectral coverage in orbital phase. The basis for this 
limitation to the range of the CfA spectra is the experience 
that, in a very eccentric orbit, the coverage of the total velocity 
range is more important than the extension of the spectral range 
(nonetheless, a significantly broader spectral range would help if 
it also included a minimum of spectra near periastron passage). 
The orbital period was fixed to the accurate value determined 
by TSL97 (Table [3]). The Elodie and Hermes spectra were 
oversampled to the velocity resolution of the CfA spectra. The 
code minimizes the squared differences between the observed 
composite spectra and the model spectra computed with respect 
to a chosen set of orbital elements. The objective function, x^i 
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Fig. 1. Uncertainty of the radial velocity information content as 
a function of wavelength, assuming S/N=100. The standard de- 
viation, Srv, was computed over a 60 A wide interval, shifted 
with a step of 10 A. The primary and secondary are represented, 
respectively by a full and a dashed line. Each wavelength region 
(Table[TJ is represented by a small horizontal line. 

can thus be written as follows: 

X^^Ya ^W[^°''^.'0') - (.^aIaU + Sa,,) + €Bh{j + 6B,,)f, (1) 
J 

where lohs.tiJ) represents a set of normalized composite spectra 
at orbital phase t and pixel j with weight W{t), IaU) and Isij) 
are the two component spectra, 5^ , and 6B,r are the component 
Doppler shifts. The component's respective (time-independent) 
light contributions are indicated by and I^b- Since a com- 
plex structure with secondary minima was revealed in the 
space (due to the broad lines of the secondary component), a 
search for the global minimum was executed for a grid of Kg- 
values (Fig. |2]l. The final spectroscopic orbit computation was 
performed with Kb - 43.4 km s"' determined from this grid 
search (cf. Table [3] upper panel). We adopted an uncertainty of 
the order of 0.5 km s ' on this value, corresponding to treating 
as equivalent all the solutions with 298840 < x^ < 298850, as 
indicated by Fig. |2] Unfortunately, Fourier disentangling does 
not (yet) include any error estimation of the orbital elements. In 
step (b), the component spectra were reconstructed applying the 
separation algorithm with this spectroscopic orbit in other wave- 
length regions (Table [T]i. 

In order to verify whether the numerical problem is well- 
conditioned, the condition numbers C,„ (i.e. the ratio of the 
largest to the smallest eigenvalue in the covariance matrix) for 
the set of equations for each Fourier mode m was computed (see 
Fig. [3]). log C is a rough measure of the loss of precision on 
the Fourier amplitudes in the component spectra, expressed in 
units of number of digits. One digit is lost, relative to the bot- 
tom value of log C, in the m - \ mode and this provoked a 
low-level sinusoidal undulation in the reconstructed spectra (see 
[Hensberge et al. 2008 for more explanation). Here, mode is 
completely undetermined (log C = oo) due to the lack of eclipses 
(dilution of spectral lines not varying with time). The low-level 
sinusoidal undulation in the reconstructed spectra due to the 
lack of cons traint on the luminosity ratio was removed self- 
consistently (llliiic et al . 2004 and Hensberge et al. 2008) This 
corresponds to replacing the best purely mathematical solution 
by the most acceptable solution from the viewpoint of physics 
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Fig. 2. Distribution of the x^ of the solutions derived by spectra 
disentangling as a function of the radial velocity semi-amplitude 
of the secondary. 




200 300 
Fourier Mode 



Fig. 3. Logarithm of the condition number (log C) versus the 
Fourier mode computed for 9~ Tau using the radial velocities 
derived with FDBinary. 



(flat continuum, conserving observed line blocking) at the ex- 
pense of a completely insignificant increase in the ;t'^-value. 

The indeterminacy of the m - mode implies that the light 
ratio between the two components has to be estimated from ex- 
ternal information. Therefore, it was determined for each spec- 
tral region using the differential magnitude (Am) measurements 
reported by AM06 (cf. their Table 3). A linear dependence of Am 
on A represents these measurements within their uncertainties. A 
weighted fit gives: 



Am = 1.131 +0.020 (/I- 6361)/1000. 
+ 0.012 + 0.011 



(2) 



The zero-point at the weighted average wavelength of 
6361 A was chosen such that the constant and the gradient are 
uncorrected. Eq. ^ fixes the light ratio (see black line in Fig. 
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Fig. 4. Line blocking of 0^ Tau AB for each of the 7 spectral 
regions described in Table[T] See Sect.|3] 
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Fig. 5. Determination of the light contributions of 6^ Tau AB, 
using the disentangled component spectra. The grey vertical bars 
represent the acceptable values for ts/tA, for each spectral re- 
gion (Table [TJ. The {b/^a predicted by AM06 are also shown 
(i.e. black diamonds with eiTorbars), as well as Eq. (|2]l (straight 
black line). The grey line represents the slope of the flux ratio 
curve of two synthetic spectra with a difiference of temperature 
corresponding to +200 K, in the sense (component B - compo- 
nent A). See Sects. l3]andl4H 



|5]l, used to renormalize the component spectra, with the central 
wavelength A - A^ . of each region /' (cf. Table [TJ substituted in 
the equation. 

In a first step, we fixed the value of the luminosity ratio in 
order to verify whether the disentangled component spectra are 
consistent with the assumption of equa l line blocking, since both 
components have very similar colours (iPeterson et al.l[T993h and 
must have the same overall chemical composition. Fig. |4] com- 
pares the line blocking for each spectral region in both compo- 
nents. The result is sensitive to sub-percent changes in the con- 
tinuum levels, in other words, to the solution of the Fourier mode 
m = 0. Physical constraints on pseudo-continuum data points de- 
liver limits to the coupled continuum placement (represented by 
grey and black symbols). The diagonal in the figure represents 
the locus of equal line blocking. As seen from this figure, the 
assumption of equal line blocking is nowhere (i.e. in none of the 
studied regions) in contradiction with the allowed range of the 
choice of the continuum levels. Under this assumption, the offset 
in continuum level of the primary's spectrum may be larger than 
the previously derived limits (see, for example, R6 in Fig. HI, 
but it will never exceed the level of 0.5%. Such an offset remains 
insignificant in the context of the present study. 

Therefore, in a second step, we evaluated which range of 
continuum positions was acceptable around the one for assumed 
equal line blocking, and, in turn, estimated the range of possible 
monochromatic luminosity ratios from the disentangled compo- 
nent spectra. Fig. |5] shows that our spectroscopic estimates are 
in concordance with the interferometric ones, though are less 
accurate than the latter. The grey intervals represent the possi- 
ble values for ^b/ ^^ derived from the renormalization procedure 
dlliiic et al.ll2004l) used to correct the disentangled component 
spectra into physically meaningful spectra. 



4. Component Spectra 

4.1. Comparison with a reference stellar spectrum 

We compared the reconstructed component spectra to the ob- 
served spectrum of the 6 Scuti star HD 2628 (Vsin / ~20 km s"' ) 
which is classified as A7 III, as is also the case of the binary 
system 9^ Tau AB (AM06). For visual concordance, the refer- 
ence spectrum was shifted by 39.5 km s ' and broadened by 66 
km s"' and 118 km s"' in order to match the spectra of compo- 
nent A and B respectively. 

Fig. |6] shows how well the reconstructed component spec- 
tra agree with the reference star spectrum, especially for the 
spectrum of component A which matches better the luminosity 
class of HD 2628. The spectra of component B may, at the one 
per cent level, still be affected by low-amplitude wiggles corre- 
sponding to larger uncertainties in the solution of the m = 2 or 
even slightly higher Fourier modes. The position and shape of 
the lines coincide well, showing the power of the spectra dis- 
entangling technique in reconstructing component spectra with- 
out any a priori assumptions about their spectral features. Note, 
however, that the shape of the Balmer Hy6-line (Fig. |6] panels 
R2) indicates a slight temperature difference between the two 
components of &^ Tau and the reference star (HD 2628 is about 
400 K cooler than 6^ Tau AB). 

4.2. Comparison with synthetic spectra 

In order to determine Vsin / and the system velocity y, we also 
compared the reconstructed spectra with synthetic spectra in the 
regions R3 to R5 of Table [T]and averaged the results. 

Independently fitting the Ha and HyS-lines of each compo- 
nent spectrum, an average effective temperature over these two 
H-lines was determined: 7800 + 170 K for component A and 
B. We also derived two limits for the temperature difference. 
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Fig. 6. Comparison between the reconstructed component spectra of 0^ Tau AB (grey) and the spectrum of HD 2628 (black), for 
the regions listed in Table[T] The latter spectrum was broadened by 66 and 118 km s respectively for components A and B, and 
shifted by 39.5 km s"'. 
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ATeif3-A, by computing synthetic spectra with slightly different 
temperatures and by imposing that their flux ratio curve should 
pass through all interferometric measurements (cf. black error 
bars in Fig.|5]). Note that this curve is more sensitive in the bluest 
part of the wavelength range. These limits suggest a temperature 
difference, ATeff^-A, ranging from +100 to +500 K (best fit at 
+200 K, cf. grey line), which is also consistent with our derived 
values of the component temperatures, considering their uncer- 
tainty. 

A small difference in effective temperature is also confirmed 
by the interferome tric measurement o f the colour difference of 
~ -0.006 mag (Pete rson et alj [1993 ). For both components, 
log g was estimated from the averaged effective temperatures 
and the components luminosities (see Sect. l5.2l i. Table|2]summa- 
rizes the atmospheric parameters of 0^ Tau using the synthetic 
spectra and compares them with the parameters adopted (except 
y^sin / and y, which were determined by TSL97). 

The Stromgren photometric d ata provided by the Ge neral 
Catalogue of Photometric Data (iMermilUod et al.l Il997h for 

Tau : V = 3.41, b-y = 0.100, oti= 0.197, ci= 1.012, u-b 
= 1.606, and fi - 2.831, confirm these effective temperatures. 
Indeed, a combined effective temperature Teff.AB = 7928 K re- 
sults from an updat e d ver sion of the standard calibration by 
iMoon & Dworetskvl (Il985h (Napiwotski, private communica- 
tion). Fig. [T] shows the disentangled component spectra and the 

Table 2. Atmospheric parameters of 6^ Tau A and B. 



Atmospheric 


Synthetic 


TSL97 


parameters 


spectra 




VAsin i (km s"') 


68.4 ± 1.5 


70^ 


Vgsin ; (km s"') 


113 + 6 


110±4 


Tcff.A (K) 


7800 ± 170' 


8250- 


Teff,B (K) 


7800 ± 170' 




log gA 


3.6 ±0.1 


4.0^ 


log gB 


3.9 ±0.1 


4.52 


y (km s~') 


39.3 ± 0.9 


39.5 ± 0.2 



Notes. The atmospheric parameters and their standard deviations were 
computed from the component spectra using sy nthetic spectra (sec ond 
column), in this work. The values determined bv lTorres et alJ ( fT997h are 
also listed for comparison (third column). 

*'' average of two values obtained using Htt and Hy6-lines; adopted 
by TSL97. 



synthetic spectra computed for the solar metallicity and for two 
possible values of the Hyades metallicity proposed by LEOl 
([Fe/H]=+0.14 and [Fe/H]=+0.19). The observed line strengths 
suggest that higher metallicities might be better than solar ones, 
as expected for two stars belonging to the Hyades cluster, but this 
claim should be confirmed by a detailed abundance analysis. 

5. Combined Orbital Analysis 

Accurate orbital solutions of stellar systems are obtained when 
different techniques are used together to provide input data to a 
single least-squares analysis, thereby leading to the simultane- 
ous determination of all the parameters, also including an inde- 
pendent determination of the distance to the system. To this aim, 
a combined interferom etric-spectroscopic analysis of 6^ Tau was 
performed dLampens^et al. 2009). Unfortunately, to our knowl- 
edge, no software is available that directly combines the orbital 



solution of spectral disentangling with the astrometry. Thus, ra- 
dial velocities must be specified, with uncertainties. These un- 
certainties will define the relative weights of the radial veloci- 
ties with respect to those of the relative positions (interferome- 
try). We used 2 x 127 radial velocities obtained from the Elodie, 
Hermes and CfA spectra and 34 best-fit angular separations p 
and position angles (Sect.|2|. 

5.1. Input radial velocities and their uncertainties 

There exist two options for the input radial velocities: (a) we 
may derive radial velocities applying the cross-correlation tech- 
nique using the reconstructed component spectra as templates, or 
(b) we may compute radial velocities from the spectroscopic or- 
bit derived with FDBinary. The first delivers "observed" (noisy) 
measurements, but the orbit derived from them will not ex- 
actly correspond to the orbit derived with FDBinary. The orbit 
that would be derived from the radial velocities resulting from 
the cross-correlation is of slightly inferior quality than the self- 
consistent orbit obtained from FDBinary. Hence, option (b) was 
preferred. Nevertheless, we checked that choosing option (a) is 
consistent with the combined astrometric-spectroscopic solu- 
tion. We also evaluated the radial velocity information content 
in the spectral region, assuming that random noise largely domi- 
nates the error budget in order to derive associated uncertainties 
(see e.g. Fig. [TJ. In this sense, these estimates are lower limits 
to the true uncertainties. The principles for a single-star spec- 
trum are elaborated in se veral papers, e.g. lVerschuer en & DavidI 
(|i999) and for a binary in lHensberge et al.U 2000). The latter un- 
certainties essentially contain a correction factor (relative to the 
single-star case) for the correlation at a given orbital phase be- 
tween the spectral gradients in the Doppler-shifted component 
spectra, and a multiplicative factor inversely pr oportional to th e 
contribution of the component to the total light dVaz et alJ2010l) . 
Typical uncertainties thus derived are given by the following me- 
dian values: cr(RVA) - 0.15 and cr(RVB) - 0.68 km s"^ (Elodie 
spectra); cr(RVA) = 0.10 and cr(RVB) = 0.54 km s"' (Hermes 
spectra); cr(RVA) = 0.65 and cr(RVB) = 3.10 km s"' (CfA spec- 
tra). In addition, we also made sure that these uncertainties in 
radial velocity are compatible with the scatter that one would 
obtain by (re)computing the radial velocities from the range of 
acceptable orbits derived with FDBinary which correspond to 
equivalent solutions in the plane (see Fig. [2]i. Table [4| lists 
the uncertainties on the component radial velocities derived with 
FDBinary. 

5.2. Final orbital elements 

The combined solution was computed using the code vbsb2, 
which performs a global exploration of the par ameter space fol - 
lowed by a local least-squares minimization (|Pourbaixlll998h . 
The essence of the method resides in a simultaneous adjustment 
of various data types performed in two steps: (a) the minimum of 
the objective function is globall y searched following th e princi- 
ples of simulated annealing (SA. Metropolis et al.lll953l) and (b) 
the best estimate of a large number of trials provides the starting 
point of a local least-squares minimization (using the procedure 
of Broyden-Fletcher-Goldstrab-Shanno). 

Table [3] (bottom panel) lists the orbital elements and their 
standard deviations obtained after 200 runs with S A, exploring a 
small interval around the period provided by TSL97. The orbital 
elements were slightly updated in the combined astrometric- 
spectroscopic solution (cf. Table [3), without significant effect 
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Fig. 7. Intrinsic spectra of 9^ Tau A and B (in black) compared to synthetic spectra computed for the Hyades metallicity 
([Fe/H]=+0.14 in red; [Fe/H]=+0.19 in green; LEOl), and the solar metallicity (in blue). The luminosity ratio is fixed to 0.35. 
The Teff, Vsin / and log g values were fixed to those ones described in the Table |2] 



on the derived component spectra. Interestingly, our determi- 
nations of the mass ratio (0.754 + 0.001) and Kg are in close 
agr eement with the first estimates derived by Pet ersonI (Il99lh 
and Pete rson et al ] (ll993l) (see Table 1 in AM06 for the mass ra- 
tio estimate using Peterson's corrected radial velocities). Table|4] 
lists the radial velocities corresponding to the orbit derived with 
FDBinary, their uncertainties (cf. Sect. 15. lb and the differences in 
the sense (Combined solution - FDBinary orbit). The combined 
astrometric-spectroscopic orbit predicts radial velocities that 
differ less than 0.1 (for comp. A) or 0.14 km s"' (for comp. B) 
from the pure spectroscopic orbit derived with FDBinary over 
the whole time interval covered by the observations and over all 
orbital phases, except for a narrow phase interval around peri- 
astron where the differences grow to 0.7 and 0.9 km s"' respec- 
tively. This corresponds to a change of ~ 1 % in the mass of the 
components. The visual (interferometric) orbit recently derived 
by AM06 comprises 7 orbital elements of which one, namely 
the orbital period, was adopted from TSL97. The remaining six 
elements (T, a, e, i, f2, w) are in very good agreement: the differ- 
ence between AM06 and the new combined solution amounts to 
0.4cr for T, 1 .4-cr for a, 2.1 -cr for e, 1 A-cr for i, 0.4-cr for oj, and 
1 .2cr for the node Q.. Note that the latter element is modified by 
180° with respect to our solution as we follow the spectroscopic 
convention (while TSL97 followed iPan et al hereafter 
PS92). The marginally significant change in the value of the 
eccentricity places it now in-between the discordant values of 
TSL97 and AM06. As an additional check, we verified that the 
residuals in both position angle and angular separation show no 



systematic offset and the standard deviations, cr(o-c)g - 1 -2° and 
C(o-c)p = 0.47 mas, agree with the published error bars on the 
measurements. 

In conclusion, the astrometric part of the combined solution 
is in excellent agreement with the pure astrometric solution com- 
puted by AM06, except (marginally) for the eccentricity. Our 
value is also constrained by the spectroscopic data, and lies in- 
between the value derived by AM06 and that of TSL97 (see 
Table [3]l. Also note that the uncertainties of the visual orbital 
elements are very similar to those published by AM06. 

The new orbital parallax has a relative error, a-„ln, of 0.5%, 
which is several times more precise than the one of TSL97 
(ji - 21.22 mas, o-„ln - 3.6%) and the ones meas ured b y the 
Hipparcos astrometric satellite (Perrvman & ESAI 119971 : n = 
21.89 mas, o-„ln - 3.8%; ivan Leeuwen ,2007; tt = 21.69 mas, 
o-jt/ti - 2.1%). Though the value of the new parallax is slightly 
smaller, it remains compatible with these formerly derived val- 
ues, within the larger uncertainties of these previous determina- 
tions. However, the new value is in d i sagree ment with the secular 
parallax derived bv lde Bruiine et all (1200 lb (see Table[3](cr;f/;7r = 
1.6%) and Sect.lOT). 

The combined orbital solution is graphically illustrated 
by Figs. [8] (astrometric observations and combined solution) 
and |9] (spectroscopic observations and combined solution). The 
astrometric-spectroscopic orbital solution previously derived by 
TSL97 is shown for comparison. We remark that these authors 
adopted some visual orbital elements of the preliminary orbit 
published by PS92 in their combined analysis. The most con- 
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Table 3. Orbital solutions of 9^ Tau A and B from various tech- 
niques. 



Table 4. 6^ Tau AB radial velocities. 





Spectroscopic Orbital Solution (FDBinary) 


HJD + 
2400000 


Orbital 
phase 


RVa 

(km s-') 


(C - 
FDB) 


RVfi 
(km s-') 


(C - Set 
FDB) 


Parameter 




This work 


47844.67396 


0.68250 


+3.217±1.208 


+0.034 


-4.264±5.456 


-0.038 1 








47942.61109 


0.37842 


-9.362±1.130 


+0.022 


+ 12.409±5.235 


-0.021 1 


P (days) 




140.72816' 


47957.58553 


0.48483 


-5.451 ±0.747 


+0.035 


+7.225±3.442 


-0.038 1 


T 




1993.0751 


48142.89616 


0.80161 


+ 11.265±0.773 


+0.001 


-14.933±3.593 


+0.005 1 


e 




0.734662 


48145.81918 


0.82238 


+ 13.202±0.701 


-0.012 


-17.498±3.262 


+0.023 1 


(^B (°) 




235.88 


48158.82114 


0.91477 


+26.806±0.697 


-0.185 


-35.530±3.274 


+0.251 1 


(km s" 


) 


32.74 















43.40" (± 0.5, cf. text) 



Interferometric Orbital Solution (AM06) 



P (days) 


140.72816' 




T 


1989.60701 + 0.00003 




e 


0.73725 + 0.00036 




a (mas) 


18.796 ± 0.056 




i(°) 


47.61 ± 0.09 




an 


173.73 ± 0.07 




<^A (°) 


55.40 ± 0.06 




T!sec (mas) 


22.50^^ ±0.36 




Astrometric-Spectroscopic Orbital Solution 


Parameter 


This work 


TSL97 


P (days) 


140.7302 ± 0.0002 


140.72816 ±0.00093 


T 


1995.00144 + 0.00002 


1993.0752 ± 0.0008 


e 


0.7360 ± 0.0003 


0.7266 ± 0.0049 


a (mas) 


18.91 ± 0.06 


18.60* ± 0.20 


i(°) 


47.8 ±0.1 


46.2-' ± 1.0 


no 


353.82 + 0.09 


171.2'' ± 1.8 


<^B (°) 


235.41 ± 0.08 


236.4 ± 1.1 


^ Ma 


0.754 ± 0.001 


0.873 ± 0.048 


T^orb (mas) 


20.90 ±0.10 


21.22 ± 0.76 


Ka (km s-') 


32.95 ± 0.04 


33.18 ±0.49 


Kb (km s"') 


43.68 ± 0.14 


38^ ±2 


Ma (Mo) 


2.86 ± 0.03 


2.42 ± 0.30 


Mb (Mo) 


2.16 ±0.02 


2.11 ±0.17 


rss (all) 


1.23e+3 


8.18e+03 


System mass (Mq) 


5.02 ± 0.09 


4.54 ±0.51 


Time span (yr) 


20.0 


6.3 



Notes. The standard deviations including orbital parallax and masses 
are shown in the bottom panel. 

<" Fixed (adopted frorn TSL97) Fixed (see Sect.^ Adopted from 
Ide Bruiine et alj d200lh (not an orbital element) Adopted from PS92 
Fixed (after a search in 2-D space). 



spicuous difference between both studies (ours - TSL97) is the 
larger radial velocity amplitude of component B, while there is 
excellent agreement for the radial velocity curve of component 
A over the entire phase range. Note that, even though a homo- 
geneous coverage in orbital phase was achieved, we still do not 
have a perfect coverage in radial velocity in the critical part of 
the orbit where the radial velocity is changing fastest (only 13% 
of the spectra have a radial velocity larger than 24.8 km s"' for 
component A). 

At this stage, it is relevant to compare and explain the un- 
certainties quoted in Table [3] Here, we list the formal errors of 
a (local) least-squares minimization process. These show a sig- 
nificant improvement of (at least) a full order of magnitude over 
the errors listed by TSL97. However, TSL97 did not compute 
their solution using full least-squares analysis (since the inter- 
ferometric data were not available to them). In addition, some 



Notes. The radial velocities (RV) were derived relatively to the center 
of mass (y = 39.3 km s"') with the code FDBinary. (C-FDB) refers to 
the differences in the sense Combined astrometric-spectroscopic solu- 
tion - FDBinary orbit. The uncertainties were computed as explained in 
Sect. 5.1. The Heliocentric Julian Day (HJD) and orbital phase corre- 
sponding to each radial velocity are also included. Set 1: CfA data; Set 
2: Elodie data; Set 3: Hermes data. A complete version of this table is 
available on-line. 



of their uncertainties were derived by way of Monte Carlo simu- 
lations assuming Gaussian errors on the various parameters. As 
akeady mentioned, with regard to the astrometric part of the so- 
lution, TSL97 adopted the visual orbital elements and the coiTe- 
sponding eiTors from PS92. Compared to this and other previous 
analyses performed with (subsets of) the Mark III data, the as- 
trometric orbit derived by AM06 shows a major improvement in 
quality: from their Table 3, a factor of about 10 was gained in the 
accuracy of the visual orbital elements over that of the previous 
studies. Our work merely confirms the very high quality of their 
orbital solution. 

On the other hand, with regard to the uncertainties in the 
spectroscopic part of the combined solution, we notice that the 
higher quality (higher S/N) of the 44 Elodie and 1 3 Hermes spec- 
tra has led to the extraction of 5-10 times more precise radial 
velocities: the quality of the radial velocities is expected to be 5 
times better if the ratio of their typical uncertainties is considered 
(e.g. Elodie versus CfA, see Sect. 15. lb . and 10 times better if the 
ratio of the typical rms residual is considered (TSL97 derived a 
rms residual of 1.7 km s^'in RVa whereas we have rms residuals 
of 0.15 and 0.20 km s Hn RV^i and RVg, respectively, over the 
entire data set). The gain might be somewhat larger considering 
the larger amount of spectra (127 instead of 70 spectra), as well 
as the robustness of the applied disentangling technique (the ra- 
dial velocities in the observed spectra are bound by a Keplerian 
orbit which probably lowers the bias in the different parame- 
ters) and the improved astrometric orbit. But the somewhat less 
favourable distribution over the orbital velocity range of the new 
spectra in comparison with the CfA ones, might have limited 
the gain achieved. Hence, an overall improvement of a factor of 
(at least) 5 can be expected compared to TSL97. As a matter of 
fact. Table |3] suggests a gain with a factor slightly larger than 10 
in the radial-velocity amplitudes. 

This difference can be understood as these formal errors re- 
main underestimates of the true errors. For example, the error 
mentioned for K^i is + 0.04 km s"\ which is about 3 times better 
than the median error in radial velocity extracted for that compo- 
nent (see Sect. 15.1b . This seems reasonable. However, the error 
mentioned for Kg is + 0. 14 km s"' , which is about 5 times better 
than the median error in radial velocity extracted for that compo- 
nent (see Sect. 15.1b . An uncertainty 2 (perhaps 3) times as large 
for Kb might be closer to the true uncertainty value, in particu- 
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Fig. 8. Combined orbital solution plotted with the astrometric 
data fi'om .Armstrong et al.. (.2006.) . 

lar if we also consider the shape of the minimum in Fig. |2] 
Indeed, the formal error based on an increase by a unit of is 

0. 15 km/s (assuming a perfect parabolic shape near minimum), 

1. e. similar to the uncertainty listed in Table [3] However, such 
an estimate corresponds to a strict lower limit of the uncertainty 
since it assumes random errors. With resampled data and the am- 
biguity in tracing of the continuum levels, uncertainties on sub- 
sequent pixels get correlated, which unavoidably contributes to 
introducing some bias in one spectrum with respect to the others. 
To stay on the safe side, we will adopt from hereon an error twice 
as large for K^, i.e. + 0.28 km s"' (instead of + 0.14 km s"'). 
As a consequence of the laws of error propagation, we will also 
consider a larger error contribution of the semi-axis major ub 
(the semi-axis major of component B with respect to the centre 
of mass expressed in km s"') and on M^i by a factor of 2. The 
former leads to an error budget on the orbital parallax increased 
by a factor of 1 .4 (determined by the ratio of the apparent and 
the true semi-axis major). In Table|5] we summarize all redeter- 
mined component properties including the increased error bud- 
get and we compare the new values with those of, respectively, 
AM06 and TSL97. With respect to AM06, the present total mass 
is 20% larger, most of which is due to the difference in parallax 
value. With respect to TSL97, the total mass is 10% larger, with 
half of the increase due to the updated parallax and the other half 
due to the increased sum of the radial velocity amplitudes. 

Adopting the new parallax of Tableland Amy = 1.11 + 0.01 
mag (derived from Eq. (I2]i) we obtain the following component 
absolute magnitudes: Mv^i = 0.33 + 0.03 mag and Mv^ - 1.44 
± 0.03 mag. This determination is more accurate than before 
(thanks to the improved parallax determination) and remains in 
good agreement with the one provided by TSL97. However, it 
differs from that by AM06 at the 3cr-level. 

6. Comparison with evolutionary models 

6.7. Previous Status 

Various attempts to compare the components' locations in the 
Hertzsprung-Russell (H-R) diagram with evolutionary models 
have akeady been done, though not always with great suc- 
cess. The importance of finding the correct evolutionary mod- 
els may also have implications for the entire Hyades clus- 
ter. iLastennet et al.l (fT999i) tested 3 independent sets of stellar 




0.4 0.6 

Orbital phase 

Fig. 9. Combined orbital solution (full lines) plotted with com- 
ponent radial velocities (symbols "o" for component A and "□" 
for component B) derived with FDBinary. Combined orbital 
solution derived by TSL97 (dashed lines). The different data 
sets obtained with the Elodie (Observatoire de Haute-Provence, 
OHP), Hermes (Roque de Los Muchachos Observatory, LPA) 
and CfA spectrographs are represented by the different shades 
of grey, see legend. 




Fig. 10. Location of 0^ Tau AB in the Hertzsprung-Russell di- 
agram. The curves represent evolutionary tracks for the masses 
corresponding to (f Tau AB, the metallicity of the Hyades (fol- 
lowing LEOl) and three levels of overshooting, orov, (different 
black curves as indicated in the legend) and an isochrone of 
age of 650 Myr (red curv e), all computed using the code Cesam 
(iMorel & Lebretonl2008l) . 



evol utionary tracks: those of the Geneva and Padova groups 
(cf. Lastennet et al. (Il999h for the references) and those of 
IClaret & Gimene J ( 1992h . Thev found some discrepant results 



from isochrone fits for 3 visual binaries using the Hipparcos 
parallaxes (including 6^ Tau) and were unable to find a sin- 
gle isochrone passing through the locations of both com- 
ponents of V818 Tau (van Bueren 22). However, LEOl re- 
marked that the models previously used (taken from the litera- 
ture) did not match the cluster's chemical composition. Using 
the location of V818 Tau, these authors were able to ex- 
clude the isochrone at 625 Myr with the solar-scaled abun- 
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Table 5. Fundamental properties of 0^ Tau A and B components. 



Component A 


This work 


AM06 


TSL97 


Teff,A (K) 


7800 ± 170 


— 


8250' 


V^sin i (km s"^) 


Do. 4 + l.J 




7ni 
/u 




1 . / / ± U.UJ 


1 7A a- n 


1 70 . n 07 
i . /o + u.u / 


IVlyi (Mq) 


9 -t- n OA 




9 zi9 -t- n "^n 


Mv/ (mag) 


U.jj ± U.Uj 


U.45 + U.Uj 


U.J / ± U.UO 


log g (cm s j 


i.b ± U. 1 




4.U 


Component B 


This work 


AM06 


TSL97 


Teff.B (K) 


7800 ± 170 




8250' 


V j^Ji-Vl. L y J\_lll J / 


113 + 6 




110 ± 4 


loe(L/Ln) 


1.32* + 0.04 


1.29 + 0.06 


1.34 ± 0.07 


Mh (Mn) 


2.16 ± 0.02 


1.87 + 0.11 


2.11 ± 0.17 


Mv (mag) 


1.44 ± 0.03 


1.61 +0.06 


1.47 ± 0.08 


loe e (cm s~^) 


3.9 ± 0.1 




4.5' 




This work 


AM06 


TSL97 


'1 Ma 


0.754 ± 0.002 


0.873^ 


0.873 ± 0.048 


Am55o„,„ (mag) 


1.11 ±0.01'' 


1.12 ±0.03 


1.10 ±0.01^ 


;T„rt (mas) 


20.90 ± 0.14 


22.30 ± 0.36^ 


21.22 ± 0.76 


(%) 


0.7 


1.6 


3.6 


Tot. mass (Mq) 


5.02 ± 0.12 


4.03 ± 0.20 


4.53 ±0.51 


IT'S,™ (%) 


2 


5 


11 


Method used 


VB-SB orbit VB orbit + ;r,,,. 


, VB-SB orbit 



Notes. The fundamental p roperties of Tau A a nd B derived in this 
wo rk (second column ), by I Armstrong et al] i2006l) (third column) and 
by [Torres et"al] ( Il997h (fourth column) are described here. The listed 
uncertainties may be slightly different from the formal uncertainties in 
Table[3]as they t ake into accoun t the increased error budget of Kg. 
(" adopted; <2' de Bruijne et al. (2001); <3) from iTorres et al] jl997l) ; 
from Eq. Q from .Peterson et al.. H993h Mb„,,o = 4.75 mag 



dance for He (Y=0.28), and derived the improved value Y - 
0.255 + 0.009. They furthermore concluded that an isochrone 
with ([Fe/H]=+0.19, Y=0.27) matched best the observed main- 
sequence of the Hyades, whereas a chemical composition with 
both ([Fe/H]=-H0.19, Y=0.27) and ([Fe/H]=-h0.14, Y=0.26) fit- 
ted well the locations of the 5 considered binaries (including 
(f- Tau). In the turnoff region, rotational effects complicate the 
interpretation but the overall agreement is good (see Fig. 12 from 
LEOl). However, in order to discriminate between the different 
models (e.g. with or without overshooting), the accuracy on the 
component masses, also of Tau, was still insufficient. 

AM06 adopted the d ynamical parallaxes b ased on the 
Hipparcos proper motions jde Bruijne et al.ll200 1') to revise the 
orbital solution of 0^ Tau (we recall that this parallax is incon- 
sistent with our determination). These authors derived a total 
mass of 4.03 + 0.20 Mq. Using the spectroscopic mass ratio of 
TSL97, they computed the component masses and luminosities. 
They were however unsuccessful at modelling the component 
locations in the mass-luminosity and the Hertzsprung-Russell 
diagrams, even though they corrected the luminosities and the 
colours for the rotational effects assuming 3 different inclina- 
tion angles: in each case, the masses appeared too small for the 
observed luminosities. Therefore, they claimed that the compo- 
nent properties of (P- Tau did not match the cuiTent evolutionary 
tracks. However, the He abundance used in their model fits is not 
appropriate for the Hyades. 



6.2. Current Status 

The revised component properties (cf. Table|5]) and, maybe more 
importantly, the higher accuracy with which they were obtained, 
make it worthwhile to review the location of the components of 
6*^ Tau in the evolutionary diagram. An isochrone model with the 
age and the chemical composition of the Hyades should be fitting 
the parameter box of both stars within the quoted uncertainties. 
We may also expect that sharper constraints might be put on the 
age and the abundance determination of the cluster, given the 
fact that the accuracies are (expected to be) of the order of 2% 
on the component masses and of the order of 4% on the compo- 
nent luminosities (i.e. an improvement with a factor of 6-7 with 
respect to previous determinations). This was already the case 
with V818 Tau (LEOl), an eclipsing-spectr oscopic binary whose 
component masses are known within 1% dPeterson & Solenskvl 
119881) . 

In Fig. (To] we plotted the locations of 6^ Tau A and B in 
the H-R diagram using the temperatures and the masses from 
Table |5] We compare these locations with theoretical evolution- 
ary tracks adopting the Hyades composition derived by LEOl 
([Fe/H] = -hO.14 and Y = 0.26). Three different values of over- 
shooting were considered: zero overshoot (full lines), and over- 
shoot values of 0.2 and 0.4, in units of the pressure scale height 
Hp (dashed, respectively dotted lines). As can be clearly seen 
from Fig. [TO] the track with M = 2.16 Mq and zero over- 
shoot passes through the box of component B (at an age of a 
600 Myr), while the tracks with M - 2.86 Mq are way too lu- 
minous compared to the location of the box representing com- 
ponent A. Nevertheless, an isochrone model of age 650 Myr and 
with the same chemical composition (thin line) passes through 
both boxes. We conclude that component B is indeed on the 
main sequence (as previous authors did), whereas the interpreta- 
tion for component A remains enigmatic: neither convective core 
overshooting nor rapid rotation can be invoked to explain the 
observed discrepancy. Taking into account gravitational darken- 
ing due to rapid rotation would shift its observed location in the 
wrong direction (its non-rotating counterpa rt would be hotter i f 
seen equator-on and fainter if seen pole-on. [premat et al.ll2005h . 
However, while non-zero overshoot would make the star evolve 
more rapidly (the star will be younger at the TAMS position), 
turbulent diffusion due to rapid rotation would make the star 
evolve more slowly (the star would be older at the TAMS posi- 
tion, Mevnet & Maeder 2000), shifting the observed location in 
the right direction. Another possibility which seems more plau- 
sible, is a change in chemical composition - in the sense of a 
higher metallicity. Interestingly, the disentangled spectra are go- 
ing in the same sense. This, however, needs to be confirmed by 
a detailed chemical analysis of the component spectra. 

Because of its location on the tip of the turnoff region, pre- 
cisely 6^ Tau A was critical to the choice of the chemical com- 
position suitable for the Hyades cluster in LEOl's work. It might 
be necessary to look into a few other solutions again. In par- 
ticular, the higher metallicity and higher helium abundance of 
the ([Fe/H] = -h0.19 and Y = 0.27) model, which fitted well 
the remainder of main-sequence stars in the cluster, should be 
re-investigated. However, such a detailed confrontation is be- 
yond the scope of the present work. Furthermore, and for the 
first time, the availability of the disentangled component spec- 
tra opens the way to an abundance determination in both com- 
ponents as if they were single stars. The resulting disentangled 
component spectra thus provide two additional "calibrator" stars 
useful in a consistent analysis of the chemical composition of 
cluster members. 
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7. Summary 

The fainter component of the system rotates rapidly such that 
its lines are continuously blended with those of the primary dur- 
ing the whole orbit. This resulted in a broad range of mass ra- 
tios published in the literature, ranging from 0.73 to 0.873. The 
spectra disentangling technique has proven, under these difficult 
conditions, to be adequate to obtain reliable orbital parameters 
as well as the intrinsic spectra of both components. These spec- 
tra, obtained from a pure mathematical technique without any 
a priori information about the component spectral features, per- 
mitted us to derive consistent and accurate atmospheric proper- 
ties such as the component effective temperatures (Teff,A=Teff3= 
7800+170 K and the respective projected rotational velocities 
(VAsin; - 68.4 + 1.5 km s"' and VBsin/ = 113 + 6 km s"'). 
The value for the mass ratio we found is: q - 0.754 ± 0.002. 
Furthermore, the component spectra were shown to closely re- 
produce the characteristics of an observed reference spectrum 
(using the single star HD 2628). From the comparison with syn- 
thetic spectra, we found that enhanced metallicity was needed 
to model Tau A and B, as expected for two members of the 
Hyades cluster. 

The combined (astrometric-spectroscopic) analysis permit- 
ted us to determine the orbital parallax (20.90 + 0.14 mas) as 
well as the component masses {Ma- 2.86 + 0.06 M© and Mb- 
2.16 + 0.02 M0) with a high accuracy. Adopting this new par- 
allax, we obtained the component absolute magnitudes Mv^ = 
0.33 + 0.03 mag and Mv^ = 1.44 + 0.03 mag. A further im- 
provement of the accuracy on the orbital parameters might be 
gained if the interval around periastron passage were more inten- 
sively covered. Indeed, at present, 76% of the spectra cover half 
the period around apastron, whereas only 24% cover half the pe- 
riod around periastron. It might also be worthwhile to extend the 
spectral range used for the application of spectra disentangling 
by one order of magnitude (with respect to the present range of 
26 A) provided that sufficient spectra obtained near the epoch of 
periastron passage are included as well. Obviously, covering this 
crucial orbital phase with more Elodie or Hermes spectra should 
lead to a significant improvement. Ideally, the distribution of the 
spectra should be homogeneous over the entire range of Doppler 
velocities. 

From the confrontation between the observed masses and lu- 
minosities and evolutionary tracks adopting the Hyades compo- 
sition derived by LEOl, we conclude that component B is on the 
main sequence. The interpretation for component A, however, 
remains problematic. We investigated a number of possibilities 
(such as convective core overshooting, fast rotation) to explain 
the observed discrepancy but it would appear from this first anal- 
ysis that a change in chemical composition - in the sense of a 
higher metallicity - might be necessary. 

The resulting disentangled component spectra open the way 
to an abundance analysis in the components of 6^ Tau as if 
they were single stars. This will be useful in the discussion 
concerning the true chemical composition of the Hyades clus- 
ter. We intend to perform such an analysis making use of the 
semi-automatic method that we developed for A and F-type stars 
tHekker et al.,2009.) . 

Acknowledgements: K.B.V. Torres, Y. Fremat and P. Lampens 
gratefully acknowledge support in the framework of the follow- 
ing projects financed by the Belgian Federal Science Policy: the 
project "Disentangled Components of Multiple Stars as Laboratories 
of Stellar Evolution" (Supplementary Researcher 2008-2009) and the 
Action- 1 project "Pulsation, Chemical Composition and Multiplicity 



in Main-Sequence A-F Stars" (Ref MO/33/018). In addition, K.B.V. 
Torres acknowledges funding from the Brazilian agencies CNPq, 
CAPES and FAPEMIG. We thank Drs. D. Pourbaix and S. llijic for 
supplying their respective codes FDBinary and vbsb2. We further- 
more thank the Optical Infrared Co-ordination Network (OPTICON, 
a major international collaboration supported by the Research 
Infrastructures Programme of the European Commission's Sixth 
Framework Programme) for the observing runs being funded by its 
Trans-National Access Programme (Ref. 2005/042 and 2006/042) as 
well as Dr. G. Torres for providing us the spectra collected at the 
Oak Ridge Observatory. The Hermes project is a collaboration be- 
tween the KULeuven, the Universite Libre de Bruxelles and the Royal 
Observatory of Belgium with contributions from the Observatoire de 
Geneve (Switzerland) and the Thuringer Landesstemwarte Tautenburg 
(Germany). Hermes is funded by the Fund for Scientific Research of 
Flanders (FWO) under the grant G.0472.04, from the Research Council 
of K.U.Leuven under grant GST-B4443, from the Fonds National de 
la Recherche Scientifique under contracts 11SN4.4506.05 and FRFC 
2.4533.09, and financial support from Lotto (2004) assigned to the 
Royal Observatory of Belgium. The Ondi^ejov Observatory is supported 
by the project with reference AV0Z10030501. We thank P. Skoda's col- 
laborators, from Ondfejov Observatory, Czech Republic, for kindly tak- 
ing 16 spectra of G- Tau. Last but not least, we thank the referee for some 
very useful suggestions. This work made use of the Simbad astronom- 
ical data base operated by the CDS, Strasbourg (France) and the ADS 
bibliography. 



References 

Armstrong, J. T., Mozurkewich, C. A., & Hummel, D. 1992, in High-Resolution 
Imaging by Interferometry, ed. J. M. Beckers & F. Merkle. 673 

Armsti-ong, J. T., Mozurkewich, D., Hajian, A. R., et al. 2006, AJ, 131, 2643 
(AM06) 

Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 
Bertelli, G., Bressan, A., Chiosi, C, Fagotto, F, & Nasi, E. 1994, A&AS, 106, 
275 

Breger, M., Garrido, R., Huang, L., et al. 1989, A&A, 214, 209 

Breger, M., Pamyatnykh, A. A., Zima, W., et al. 2002, MNRAS, 336, 249 

Claret, A. & Gimenez, A. 1992, A&AS, 96, 255 

de Bruijne, J. H. J., Hoogerwerf, R., & de Zeeuw, R T. 2001, A&A, 367, 1 1 1 
Fremat, Y., Lampens, R, Van Cauteren, P., & Robertson, C. W. 2006, Memorie 

delta Societa Astronomica Italiana, 77, 174 
Fremat, Y, Zorec, J., Hubert, A., & Floquet, M. 2005, A&A, 440, 305 
Girai-di, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 
Hekker, S., Fremat, Y, Lampens, R, et al. 2009, MNRAS, 396, 1689 
Hensberge, H., llijic, S., & Torres, K. B. V 2008, A&A, 482, 1031 
Hensberge, H., Pavlovski, K., & Verschueren, W. 2000, A&A, 358, 553 
Hummel, C. A. & Armstrong, J. T. 1992, in ASP Conf. Ser 32: lAU Colloq. 
135: Complementary Approaches to Double and Multiple Star Research, ed. 
H. A. McAlister & W. I. Hartkopf, 552 
llijic, S., Hensberge, H., Pavlovski, K., & Freyhammer, L. M. 2004, 
in Astronomical Society of the Pacific Conference Series, Vol. 318, 
Spectroscopically and Spatially Resolving the Components of the Close 
Binary Stars, ed. R. W. Hilditch, H. Hensberge, & K. Pavlovski, 111 
Kennelly, E. J., Walker, G. A. H., Catala, C, et al. 1996, A&A, 313. 571 
Krolikowska, M. 1992, A&A, 260, 183 

Lampens, P., Torres, K., Fremat, Y, & Hensberge, H. 2009, in American Institute 
of Physics Conference Series, Vol. 1170, American Institute of Physics 
Conference Series, ed. J. A. Guzik & P. A. Bradley, 446^48 
Lastennet, E., Valls-Gabaud, D., Lejeune, T., & Oblak, E. 1999, A&A, 349, 485 
Lebreton, Y, Fernandes, J., & Lejeune, T. 2001, A&A, 374, 540 (LEOl) 
Li, Z., Zhou, A., & Yang, D. 1997, PASR 109, 217 
Mermilliod, L, Mermilliod, M., & Hauck, B. 1997, A&AS, 124, 349 
Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., & Teller, 

E. 1953, J. Chem. Phys., 21, 1087 
Meynet, G. & Maeder, A. 2000, A&A, 361, 101 
Moon, T. T. & Dworetsky, M. M. 1985, MNRAS, 217, 305 
Morel, R & Lebreton, Y. 2008. Ap&SS, 316, 61 

Pan, X., Shao, M., & Colavita, M. M. 1992, in Astronomical Society of 
the Pacific Conference Series, Vol. 32, lAU Colloq. 135: Complementary 
Approaches to Double and Multiple Star Research, ed. H. A. McAlister & 
W. I. Hartkopf, 502 (PS92) 

PeiTyman, M. A. C, Brown, A. G. A., Lebreton, Y, et al. 1998, A&A, 331, 81 



12 



Torres, K. et al.: Orbit, orbital parallax and masses of ^ Tau AB 



Ferryman, M. A. C. & ESA, eds. 1997, ESA Special Publication, Vol. 1200, 
The HIPPARCOS and TYCHO catalogues. Astrometric and photometric star 
catalogues derived from the ESA HIPPARCOS Space Astrometry Mission 
Peterson, D. M. 1991, in Astronomical Society of the Pacific Conference Series, 
Vol. 13, The Formation and Evolution of Star Clusters, ed. K. Janes, 592-597 
Peterson, D. M. & Solensky, R. 1988, ApJ, 333, 256 
Peterson, D. M., Stefanik, R. P, & Latham, D. W. 1993, AJ, 105, 2260 
Poretti, E., Buzasi, D., Laher, R., Catanzarite, J., & Conrow, T. 2002, A&A, 382, 
157 

Pourbaix, D. 1998, A&AS, 131, 377 

Press, W. H., Teukolsky, S. A., Vetterling, W. T, & Flannery, B. P 1992, 

Numerical recipes in FORTRAN. The art of scientific computing, ed. W. H. 

Press, S. A. Teukolsky, W. T. Vetterling, & B. P. Flannery 
Raskin, G., Van Winckel, H., Hensberge, H., et al. 2010, A&A, submitted 
Shao, M., Colavita, M. M., Mines, B. E., et al. 1991, in Bulletin of the American 

Astronomical Society, Vol. 23, Bulletin of the American Astronomical 

Society, 830 

Skoda, P. & Slechta, M. 2002, Publications of the Astronomical Institute of the 

Czechoslovak Academy of Sciences, 90, 22 
Torres, G., Stefanik, R. P, & Latham, D. W. 1997, ApJ, 485, 167 (TSL97) 
van Leeuwen, F 2007, A&A, 474, 653 

Vaz, L. P. v., Torres, K. B. V., Hensberge, H., et al. 2010, A&A, in preparation 

Verschueren, W. & David, M. 1999, A&AS, 136, 591 

Yildiz, M., Yakut, K., Bakr§, H., & Noels, A. 2006, MNRAS, 368, 1941 



13 



Torres, K. et al.: Orbit, orbital parallax and masses of (p- Tau AB, Online Material p 1 
Table 4. ^ Tau AB radial velocities. 
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-0.837 
-0.907 
-0.678 
-0.246 
-0.031 
+0.018 
+0.050 
+0.071 
+0.062 
+0.058 
+0.037 
+0.033 



1 50138 

1 50140 

1 50143 

1 50146 

1 50152 

1 50154 

1 53431 

1 53431 

1 53642 

1 53657 

1 53667 

1 53686 

1 53690 

1 53711 

1 53724 

1 53728 

1 53739 

1 53747, 

1 53749 

1 53759 

1 53777, 

1 53787 

1 53793 

1 53801 

1 53801 

1 53801 

1 53801 

1 53801 

1 53802 

1 53802 

1 53802 

1 53803 

1 53803 

1 53803 

1 53803 

1 53803 

1 53803 

1 53803 

1 53803 

1 53803 

1 53806 

1 53806 

1 53806 

1 53806 

1 53806 

1 53806 

1 53806 

1 53806 

1 53806 

1 53806 

1 54905 

1 54905 

1 54905 

1 54906 

1 54906 

1 54906 

1 55079 

1 55079 

1 55080 

1 55080 

1 55080 

1 55082 

1 55082 
1 



.48966 
.56345 
.47618 
47843 
.51706 
.54631 
26650 
.27540 
.66560 
52560 
.67700 
.56880 
58560 
46230 
.36150 
37360 
43640 
.35580 
39120 
35790 
.34740 
30460 
32650 
.33250 
34800 
36060 
.37240 
.38420 
30060 
.34250 
.36410 
,27240 
.28070 
.29010 
29870 
.30710 
.31560 
32400 
33260 
.34150 
29510 
30400 
.31280 
32120 
32980 
.33840 
34690 
35260 
.37650 
38500 
41312 
.42068 
.42824 
39330 
.39728 
.40149 
70896 
.71082 
.70855 
71054 
.71256 
.70400 
70576 



0.98189 
0.99662 
0.01732 
0.03865 
0.08156 
0.09598 
0.37968 
0.37975 
0.88184 
0.98743 
0.05957 
0.19381 
0.22235 
0.37070 
0.46236 
0.49087 
0.56948 
0.62575 
0.64021 
0.71103 
0.83886 
0.90962 
0.95241 
0.00930 
0.00941 
0.00950 
0.00958 
0.00966 
0.01618 
0.01647 
0.01663 
0.02308 
0.02314 
0.02321 
0.02327 
0.02333 
0.02339 
0.02345 
0.02351 
0.02357 
0.04456 
0.04462 
0.04469 
0.04475 
0.04481 
0.04487 
0.04493 
0.04497 
0.04514 
0.04520 
0.85467 
0.85472 
0.85478 
0.86163 
0.86166 
0.86169 
0.09318 
0.09319 
0.10028 
0.10030 
0.10031 
0.11446 
0.11447 



+46.222 
+36.590 
+2.621 
-13.453 
-19.179 
-19.225 
-9.296 
-9.293 
+20.691 
+44.782 
-18.027 
-16.168 
-15.083 
-9.619 
-6.272 
-5.199 
-2.069 
+0.414 
+1.099 
+4.857 
+14.995 
+25.804 
+37.826 
+14.200 
+14.010 
+13.857 
+13.713 
+13.571 
+3.621 
+3.228 
+3.028 
-4.125 
-4.179 
-4.241 
-4.297 
-4.351 
-4.406 
-4.461 
-4.516 
-4.572 
-15.442 
-15.458 
-15.474 
-15.489 
-15.506 
-15.522 
-15.537 
-15.547 
-15.590 
-15.605 
+16.850 
+16.860 
+16.870 
+17.750 
+17.760 
+17.760 
-19.240 
-19.240 
-19.180 
-19.180 
-19.170 
-18.900 
-18.900 



±1.407 
±0.605 
±0.746 
±0.677 
±0.775 
±0.630 
±0.121 
±0.121 
±0.097 
±0.257 
±0.097 
±0.097 
±0.097 
±0.145 
±0.096 
±0.072 
±0.120 
±0.240 
±0.120 
±0.096 
±0.121 
±0.123 
±0.101 
±0.145 
±0.145 
±0.169 
±0.169 
±0.145 
±0.830 
±0.265 
±0.708 
±0.145 
±0.120 
±0.120 
±0.349 
±0.145 
±0.145 
±0.120 
±0.265 
±0.651 
±0.194 
±0.194 
±0.194 
±0.218 
±0.218 
±0.266 
±0.266 
±0.691 
±0.378 
±0.266 
±0.069 
±0.062 
±0.078 
±0.081 
±0.078 
±0.078 
±0.095 
±0.095 
±0.095 
±0.095 
±0.119 
±0.095 
±0.095 



-0.466 
+0.265 
+0.520 
+0.093 
-0.058 
-0.059 
+0.042 
+0.042 
-0.001 
-0.463 
-0.070 
-0.005 
+0.005 
+0.042 
+0.055 
+0.057 
+0.063 
+0.066 
+0.067 
+0.065 
+0.038 
-0.040 
-0.221 
+0.041 
+0.031 
+0.057 
+0.047 
+0.038 
+0.055 
+0.072 
+0.075 
+0.021 
+0.038 
+0.024 
+0.015 
+0.032 
+0.024 
+0.016 
+0.031 
+0.021 
-0.056 
-0.059 
-0.055 
-0.057 
-0.060 
-0.056 
-0.058 
-0.055 
-0.059 
-0.061 
+0.017 
+0.020 
+0.020 
+0.011 
+0.018 
+0.014 
-0.058 
-0.058 
-0.059 
-0.059 
-0.050 
-0.044 
-0.044 



-61.267+8.022 
-48.499±2.979 
-3.474±3.329 
+17.831+3.169 
+25.422+3.688 
+25.482+2.996 
+12.321+0.559 
+12.318+0.559 
-27.425+0.458 
-59.357+1.464 
+23.893+0.460 
+21.431+0.457 
+19.992+0.454 
+12.751+0.674 
+8.314+0.445 
+6.891+0.333 
+2.743+0.537 
-0.549+1.044 
-1.457+0.525 
-6.438±0.437 
-19.876±0.565 
-34.202±0.584 
-50.137±0.510 
-18.821±0.676 
-18.570±0.676 
-18.367±0.789 
-18.177±0.788 
-17.987±0.676 
-4.800±3.755 
-4.278±1.196 
-4.013±3.197 
+5.467+0.662 
+5.539+0.552 
+5.621+0.552 
+5.695+1.599 
+5.768+0.662 
+5.840+0.664 
+5.912+0.554 
+5.986+1.220 
+6.061+2.997 
+20.468+0.913 
+20.489+0.913 
+20.511+0.913 
+20.532+1.028 
+20.552+1.028 
+20.573+1.256 
+20.593+1.256 
+20.608+3.262 
+20.665+1.784 
+20.684+1.256 
-22.340+0.388 
-22.350+0.348 
-22.360+0.442 
-23.530+0.456 
-23.540+0.443 
-23.540+0.443 
+25.500+0.537 
+25.500+0.537 
+25.420+0.537 
+25.420+0.537 
+25.420+0.671 
+25.050+0.537 
+25.050+0.537 



+0.622 
-0.346 
-0.681 
-0.115 
+0.086 
+0.087 
-0.048 
-0.047 
+0.008 
+0.618 
+0.101 
+0.016 
+0.002 
-0.047 
-0.064 
-0.068 
-0.076 
-0.080 
-0.081 
-0.079 
-0.044 
+0.059 
+0.298 
-0.048 
-0.035 
-0.069 
-0.056 
-0.043 
-0.066 
-0.088 
-0.092 
-0.020 
-0.043 
-0.024 
-0.013 
-0.034 
-0.024 
-0.014 
-0.033 
-0.019 
+0.083 
+0.087 
+0.082 
+0.085 
+0.088 
+0.082 
+0.085 
+0.082 
+0.087 
+0.089 
-0.021 
-0.023 
-0.020 
-0.011 
-0.016 
-0.012 
+0.083 
+0.083 
+0.085 
+0.085 
+0.085 
+0.066 
+0.066 



Notes. The radial velocities (RV) were derived relatively to the center of mass (y = 39.3 km s"') with the code FDBinary. (C-FDB) refers to 
the differences in the sense Combined astrometric-spectroscopic solution - FDBinary orbit. The uncertainties were computed as explained in 
Sect. 5.1. The Heliocentric Julian Day (HJD) and orbital phase corresponding to each radial velocity are also included. Set 1: CfA data; Set 2: 
Elodie data; Set 3: Hermes data. 



